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Abstract 

In this paper, we address the role of electron-electron interactions on the velocities of spin and charge 
transport in one-dimensional systems typified by conjugated polymers. We employ the Hubbard model 
to model electron-electron interactions. The recently developed technique of time dependent Density 
Matrix Renormalization Group (tdDMRG) is used to follow the spin and charge evolution in an initial 
wavepacket described by a hole doped in the ground state of the neutral system. We find that the charge 
and spin velocities are different in the presence of correlations and are in accordance with results from 
earlier studies; the charge and spin move together in noninteracting picture while interaction slows down 
only the spin velocity. We also note that dimerization of the chain only weakly affects these velocities. 



1 Introduction 

Low-dimensional many-body systems have always been the focus of theoretical and experimen- 
tal interest. The physics of these systems is quite different from those of three (3D) systems. For 
example, these materials show the phenomena of spin-charge separation, wherein the spin and 
charge degrees of freedom of the electron get decoupled and evolve independently of each other 
with different velocities. These materials find wide scale applications in the field of molecu- 
lar electronics (spintronics). Amongst low-dimensional materials, the ^-conjugated polymers 
have attracted a lot of interest, being potential candidates for various molecular electronics and 
spintronics applications; examples include the organic light emitting diodes (OLEDS), organic 
semiconductors, organic thin film transistors, etc. 01 SB Hi- However, spin and charge trans- 
port in these systems is still not well established because of the strong electron-electron cor- 
relations that exist in these systems. Transport in these materials is strictly a non-equilibrium 
phenomena to understand which, one needs to investigate the time evolution of strongly in- 
teracting quantum many body systems. Recently, there has been a considerable progress in 
investigation of non-equilibrium time evolution of many body systems. Analytical approaches 
like the perturbative Keldysh formalism is restricted to a few integrable models only, but in 
the case of low-dimensional systems, efficient numerically accurate techniques have been de- 
veloped and successfully applied to a variety of models. One such efficient method which has 
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gained tremendous impetus in recent years, is the time-dependent Density Matrix Renormaliza- 
tion Group technique (tdDMRG) 10,0, 0, & Oil EH]. In this paper, we use tdDMRG to address 
the effects of (1) electron-electron correlations and, (2) dimerization, on the charge (spin) trans- 
port in quasi 1 -dimensional strongly correlated polyene chains. Moreover, we also look into the 
dynamics of spin-charge separation in these systems. To address the above questions, we focus 
our attention on the real time quantum dynamics of an hole with up spin injected at site- 1 of 
polyene chains. 



2 Model Hamiltonians and Parameters 

We have modeled the ^-conjugated chains using three model Hamiltonians: (a) the Tight 
Binding (TB) Hamiltonian, known as Hiickel model to chemists lll2lll3n. (b ) the single-band 



Hubbard model QEHEl], and (c) Pariser-Parr-Pople (PPP) model ll 1 vL 1 1 8h - Amongst these, 



(2) and (3) are interacting Hamiltonians that include explicit electron correlations and (3) is a 
realistic model Hamiltonian used for describing ^-conjugated polymers. In second quantized 
formalism, these three model Hamiltonians can be expressed as given below 111911 : 

N-l 

#Hiickei = ?o £ (a] +lc a ia + h.c) (1) 
i=l,o=t,| 

N 

Hmibbard = #Hiickel + U ^ n^fl^ (2) 

i=l 

Hppp = H Hubbard + Vjj (m - zi) (n j - zj) (3) 

i>j 

Here, aj a (a !CT ) creates (annihilates) an electron at site-i of the polyene chain, to is the nearest- 
neighbour (nn) hopping integral for an undimerized chain, and h.c refers to the hermitian 
conjugate. In the case of a dimerized polyene chain, the nn hopping integral is given by, 
= fo(l — ( — S)'), where 8 is the dimerization parameter. For the present study, we have 
taken 8 = 0.07, so that the nn hoping term for long and short bonds are respectively given 
by, hong-bond = 1-07 and t s h or t-bond = 0.93, to = 1.0 for the Hubbard model, (n^) are the 
number density of upspin (downspin) electrons at site-i of the polyene chain. The Hubbard 
model is characterized by U, the Hubbard parameter, which represents on-site Coulomb repul- 
sion between two electrons of opposite spins occupying the same site of the polyene chain. For 
homogenous systems, this parameter is same for all sites. U is measured in terms of to, and the 
parameter, (U /to) characterizes the TT-electronic motion in single band systems. In our study, 
we have taken U/t = 0.0 (the Hiickel model), 2.0, 4.0, 6.0 and 10.0. In the PPP model, the Vy 
is the inter-site Coulomb repulsion between two different sites, i and j, of the polyene chain. In 
keeping with the spirit of phenomenology associated with the PPP Hamiltonian, the inter- site 
electron repulsion integrals, Vy are interpolated smoothly between U for zero separation and 
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j- for the inter-site separation tending to infinity; thus, the explicit evaluation of the repulsion 
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integrals is avoided. There are two widely used interpolation schemes used to evaluate Vu, the 



Ohno scheme Q20Q , and the Mataga-Nishimoto scheme 112 111 . In the Ohno interpolation scheme 
which we use, the inter-site electron repulsion integrals, V/ ; are given by, 



Vij = 14.397 



28.794 \ 2 2 1~ 1/2 



(4) 



The Mataga-Nishimoto formula is given by, 

2.0 



V;, 



+ ■ 



'tj 



Ui + Ui 14.397 



-l 



(5) 



The Ohno interpolation formula decays more rapidly than that of Mataga-Nishimoto scheme. 
In both the above interpolation schemes, it is assumed that is in A, while to, U, and Vij are in 
measured in ev. Zi is the chemical potential of site-i; the function of Zi is to keep the i r/, -carbon 
atom neutral when singly occupied. 



3 Time-dependent DMRG - Xiang's Algorithm 

For carrying out quantum dynamics of the an up spin hole injected at site-1 of the polyene chain, 
we first create the necessary initial state, by annihilating an up spin electron from site-1 from 
the ground state of the TT-conjugated chain. Mathematically, this amounts to the following: 

| \]/(0)) = ai t | Mfcs) (6) 



Here, | y(0)} is the desired initial state and | \\fcs) is the ground state of the neutral polyene 
chain. Using V|/(0), we numerically solve the time-dependent Schrodinger equation (TDSE) 
which is given by, 

B&m-HM)) (7) 

where H is any of the three time-independent model Hamiltonians discussed in sec. II. The 
above equation has the formal solution, 

\y(t)) = e- iHt \y(0)) (8) 



Numerically, given a small time step At, H and | y(0)), the TDSE can be solved by expanding 
the exponential function in equation (8), to various orders of (H At). The simplest of this is the 
Euler (EU) scheme as given below: 

|y(f + A0) = (l-iffAf) \y(t)}+0((HAt) 2 ) (9) 
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This equation is then repeatedly used to obtain to propagate the initial wave-packet. This 
scheme is an explicit one without requiring any matrix inversion. However, it suffers from 
two serious drawbacks: (1) it is non-unitarity, and (2) there is an instability due to the lack of 
time inversion symmetry (t — > -t). To avoid these problems, the TDSE is solved using the 



implicit Crank-Nicholson (CN) scheme 112211 in which the exponential function is approximated 
by the Caley transform 

1 -iHAt/2 . / \ \ ^ 
I + = l+iHAt/2 I V(0> +°(( HAt ) ) ( 10 ) 

The CN scheme is unitary, unconditionally stable, and accurate up to (HAt) 3 . However, this 
scheme also has a serious limitation, namely: Each time evolution step requires a matrix inver- 
sion, which for large systems and with the increase of dimensionality, requires huge memory 
and CPU time, making this method prohibitive. Hence, there has been a surge towards the de- 
velopment of explicit, stable integration schemes. The first of these is a symmetrized version 
of the EU scheme, called the second order differencing scheme (MSD2) [|23|]. This scheme is 
symmetric in time as seen below, is conditionally stable, and accurate upto (HAt) 2 . 

| y(t + At)) = -liHAt | y(0)+ I V(t-At)) + 0((HAt) 3 ) (11) 



The MSD2 scheme can be extended to higher order accuracy forms, which are collectively 



called the multistep differencing (MSD) schemes [|24ll . for example, the fourth and sixth order 



MSD (MSD4 and MSD6) which can given as below (equations (12) and,(13)): 

| + 2At) ) = I \\f(t -2At)) + AiHAt 
+ 0((HAt) 5 ) 

| \\f(t + 3At)) =| y(t - 3At)) - UHAt 



(12) 



tkIy(0>-^( lY(f+AO)+|v(f-Af)) 
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(13) 
+ 0((HAt) 7 ) 



These higher order schemes are explicit and conditionally stable, for example, MSD4 is stable 
if and only At < 0.4, while for MSD6, stability exists for At < OA. Predictor-Corrector (PC) 
techniques are another class of ordinary differential equation solvers. For our present studies, 
we have developed a PC scheme of our own, which we call the MSD4-AM4 method. In this, 
we use the explicit MSD4 (equation (12)) scheme as the predictor, and fourth order implicit 
Adams-Moultan method as the corrector (equation (14)). We found this scheme to be very ro- 
bust, and as efficient as the CN method; moreover, this PC technique is much faster and less 
memory consuming compared to the CN scheme Q25[| . 
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\\f(t + 2At)) =| \\f(t + At)) 



MAt 
24 



9 | y(t + 2A0) + 19 | y(f + At)) - 5 | V|/(f)> 



+ | Af)> +O((//A0 5 ) 



(14) 



So far we have discussed about model Hamiltonians, preparing the initial state, and time evolu- 
tion of this state by solving the TDSE numerically. For obtaining the initial and ground states of 
the polyene chains, we use tdDMRG as given by Xiang and coworkers Q26Q . However, before 
discussing this technique, we'll briefly discuss the conventional infinite system DMRG method 
[S, 5| as proposed by White and others. The basic idea of the DMRG method is to divide a given 
finite system into two parts, namely, system and surrounding, followed by retaining only the m 
most highly weighted eigenstates of the reduced density matrix of these partial "systems" Ir27ll . 
Using these reduced density matrices, one or more pure states of this total system is obtained. 
In case of the infinite system DMRG algorithm, the system size is increased in units of "two 
sites" (see fig. 1)0,0,11. 



system i surrounding 



system nl^nl ' surrounding 



system nl nl' surrounding 



RG— stepl : Increase number of 
degrees of freedom: Add 
exact site to old system 

block 



obtain |(h> 



obtain 


P 


■ 


Diagonal 


ize 


P 


1 





system 






surrounding 


nl nl' 





New Basis 


: eigenstates 


of P 




RG-step2 : 


Decrease the 


number of 


degrees of 


freedom : 


cutoff after 


m states 





Transform system block into new 
basis with only m states 



Figure 1 : Schematic sketch of infinite system DMRG algorithm along with the flowchart showing basics 
of the DMRG truncation scheme. For more details of the method, see tA\m- nl and nl' are the two new 
sites that are attached at each infinite system DMRG step. 



5 



So far, all tdDMRG schemes can be categorized into three classes: (1) Static time-dependent 
DMRG, (2) Dynamic time-dependent DMRG, and (3) Adaptive time-dependent DMRG method 
|@]. The static tdDMRG method was first introduced by Cazalilla and Marston [28], who ex- 



ploited this technique for investigating time-dependent quantum many-body effects. They stud- 
ied a time-dependent Hamiltonian, H(t) = H(0) + V(t), where V(t) represents the time-dependent 
part of the Hamiltonian. Initially, infinite system DMRG method was used for constructing a lat- 
tice of desired size keeping a substantially large number of reduced density matrix eigenstates 
(m). Time evolution of this final lattice system is then carried out using the time-dependent 
Hamiltonian, H e ff(t), which is given by, H e ff(t) = H e ff(0) + V e ff(t), where H e ff(0) is the fi- 
nal superblock Hamiltonian approximating H(0), and V e f/(t) is an approximation of V(t), and 
is built using the representations of operators in the final block bases. The basic idea of this 
method is to fix the reduced Hilbert space at its optimal value at time t = 0, and then, projecting 
all wavefunctions and operators on to it. In other words, the effective Hamiltonian which has 
been obtained by targeting the ground state of the t = Hamiltonian is capable of representing 
adequately the time-dependent states that will be reached at later times. The major disadvan- 
tage of this scheme is that it fails completely for long time evolution as there is a significant 
loss of information due to the 'final superblock truncation'. Moreover, the number of DMRG 
states, m, grows with the simulation time as they need to incorporate a constantly increasing 



number of nonequilibrium states. To overcome this, in 2003, Luo, Xiang and Wang 112611 came 
up with a targeting method, which is called the Dynamic tdDMRG or the LXW method, and 
will be utilized by us, for the present study. We will however, not discuss the Adaptive tdDMRG 
scheme. Interested readers can refer the relevant articles [|9, 10, IT, 29, 30, 31]. The algorithm, 
as implemented by us, is given in details below: 

(1) The Hamiltonian of a small, exactly diagonalizable superblock(SB) of L (= 4) sites, H[^ 4 , 
is first constructed. 

(2) The ground state, \^ s , of this 4-site SB is obtained by exact diagonalization of #££ 4 . Using 
\|/^, a desired initial state \|/q of interest is prepared. Exact time evolution of this initial state is 
then carried out from t = 0, to t = N steps , by solving the TDSE numerically, using a convenient 
integration scheme. (In our case, we use our own MSD4-AM4 scheme). At the end of the time 
evolution, a set of time-dependent wavefunctions are obtained, {v|/(//) : V G (0,A^ e/ „)}. 

(3) Using this set of time-dependent wavefunctions, the reduced density matrices for the left- 
(p/) and right-half (p r ) blocks for the next SB is build using LXW prescription jji^l . Mathemat- 
ically, 

Nsteps 

Y,®i =1 ( 15 ) 

i=0 

X>« = 1 (16) 

;=0 

Here, | is called the ith-target state and co, is its corresponding weight in the half -block 

reduced density matrix. In the original LXW method, in building of p/ and p,-, only \|/(0) 
and V i E (l,N steps ) are included. However, in our case, we have two systems at hand: 
neutral polyene, and "cationic" polyene, having +1 charge on it. We found that in the case of 
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homogeneous systems, devoid of heteroatoms, it is un-important whether we keep the "ionic" 
ground state for building of the density matrices. However, for systems with heteroatoms, 
including the "ionic" ground state is essential for building the density matrices. Furthermore, 
we have also found that by comparing tDMRG results with exact results, for small chains, eo !on 
< CGq is required. Hence, the above pair of equations get modified as, 



^steps 



p/ = TrA (0 gs j 0n I \\f gs ,ion)(y g s,ion | + £ ©i | V(*/))M f i 



;=0 



Nsteps 





^steps 






(tigs, ion 


+ E 


(0, 


= 1 




!=0 








Nsteps 






(tigs, ion 


+ E 


(0, 


= 1 




;=0 







p r = Tri [togsjon | Vgs,ion) (Vgsjon \ + £ G>i | Y(f/)) M*0 I J ! <*Wan + (0/ = 1 (18) 

(4) These reduced density matrices are then diagonalized using a dense matrix diagonalization 
routine to obtain m eigenvectors, with largest eigenvalues. These eigenvector constitute the 
Density Matrix Eigen- Vector (DMEV) basis. 

(5) Hi + \ (H r+ i) and other operators (A/+i) are then constructed in the new system block and 
transformed to the reduced DMEV basis using the transformations, Hi + \ = 0\Hi + \Ol, = 
o\Ai+\Ol- Here, Ol is a (4m x m) matrix whose columns contain m highest eigenvectors of p; 
(p r ), and A/+i is an operator in the system block (left-, or right-half block). 

(6) A new SB of size (L+2) is formed, using Hi+i, two newly added sites, and H r+ \. 

(7) The steps from (2) to (6) are repeated to iteratively increase the SB size by two sites at a 
time. 

Apart from the "ionic ground state" correction to the original LXW Algorithm, we have also 
introduced another modification, which we call the "n-slot" modification. This modification 
basically means that instead of storing all the time-dependent non-equilibrium wavefunctions, 
after every "n-th" time step, the wavefunction is stored for building the density matrix. We have 
studied "n" = 10, 100, 500 and 1000 cases. It is found that for getting correct results using 
LXW Dynamic tdDMRG technique, we need "n" > 25. The basic idea behind this method is 
that, the time-dependent wavefunctions for a SB of size L, explores the Hilbert space as much 
as possible, and transfers the information through the reduced density matrix, towards building 
Hilbert space of a SB of size, (L+2). However, this technique suffers from a major problem, 
namely, it needs large CPU times. Parallelizing this algorithm would mitigate this drawback. 
The dynamical variables that we study are charge (spin) densities at site-1, site-L of the polyene 
chain, along with charge (spin) velocities. These variables will be discussed in detail in the next 
section. 



4 Results and Discussion 

In the previous section, dynamical variables that are calculated in this paper, were mentioned. 
Here, we discuss these quantities in detail, along with our results. The charge (spin) density at 
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the F- 



site of a given polyene chain, at time x is given by: 



(n/Cc)) = (\|/(x) | (n ;t + n ;; ) 



¥(x)) 



(19) 



(5f(x)) = (x|/(x)|(n ! - t -n /; ) 



(20) 



We have calculated these two quantities at all sites of a given polyene chain. However, we 
focus our attention on (nz,(x)} ((S£(x)), "L" being the last site of the polyene chain. We have 
considered chains containing 10, 20, 30 and 40 C-atoms in the present study. Two different 
evolution times have been used namely, 33 fs (femtoseconds) and 10 fs. The dimension of the 
DMEV Basis is kept at an optimal value of 200 for the Hiickel and Hubbard chains. Site-L of 
the polyene chain at time t = has (hl) = 1 .0, with equal probability for being occupied by 
either an "up (down) spin", thereby making (5^) = 0.0. As time progresses, the injected hole 
propagates from the first to the last site, and this is represented by appearance of a minima in 
the time evolution profiles of both charge and spin densities. The time at which the 1st minima 
appears is therefore the time taken by the injected hole to reach the other of the chain. Hence, 
we focus our attention on this quantity throughout our studies. 



4.1 Dynamics in Hiickel Chains: 

Fig. 2 shows the time evolution of charge (spin) densities at the last site of polyene chains of 
different lengths, governed by the Hiickel Hamiltonian. The left- and right-plot depicts charge 
and spin density dynamics respectively. Solid curves are used for undimerized (regular) chains, 
while dashed curves, for the dimerized chains, with 8 = 0.07. In case of Hiickel chains, from 
fig. 2 it is seen that with increasing chain length, the time taken by the injected hole to reach the 
end of the chain also increases. The velocity of the hole appears to be reasonably constant for 
systems of different length. Furthermore, dimerization appears to slightly decrease the "hole- 
velocity" compared to the uniform chain. Careful examination of the time profiles of charge 
(spin) densities also reveal that they are identical, in features indicating that there is no spin- 
charge separation. 



8 



Time evolution of <n L (t)> fox Huckel chains 
|aiz = 2 Q ) r 5 = 0.07 (dotted curve) ; 8 = 0.0 (solid curve) 



Time evolution of <S 1 (t) > for Huckel chain 
(sunt = 20Q), £ = 0,07 (dnttad curve); S =0.0 (solid curve) 
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Figure 2: Time evolution of (Left) and (Sx(0) O^ght) f° r uniform (solid curve) and dimerized 

(dashed curve) Huckel chains of length 10, 20, 30 and 40 sites. 200 DMEV states per block is used. 
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4.2 Dynamics in Hubbard Chains: 



Figs. 3, 4, and 5 gives the time evolution profiles of charge and spin densities at the last site of 
Hubbard chains for different chain lengths, and for several representative values of U, namely, 
U/t = 2.0, 4.0, and 6.0, respectively. It is observed that the charge and spin dynamics are no 
longer identical as was seen in case of Hiickel chains clearly indicating spin-charge separation. 
Furthermore, as the magnitude of U increases, the extent of spin-charge decoupling also in- 
creases. In the literature, the decoupled spin and charge excitations are referred to as spinous, 



and holons 11321. 133N. For a given U/t, "velocity" of the charge excitation (holon) as well as that 
of the spin excitation (spinon) seem to be weakly dependent on the chain length. It is however 
that at all chain lengths and nonzero U/t values, "velocity" of the holon is higher than that of 
the spinon. If one examines the plots carefully, another interesting observation can be made; 
in the correlated models, dimerization plays little or no role in influencing the "velocity" of the 
injected hole. Figs. 6 and 7, show the variation of («l(x)) and (S Z L (%)) for regular and dimerized 
chains of given length, for different values of U. Solid curves are for U = 2.0, dashed curves 
represents U = 4.0, while U = 6.0 is depicted by dotted curves. It is observed that for a fixed 
chain length, increasing U, does not perceptibly affect the velocity of the holon appreciably, but 
spinon's velocity is considerably altered. This is simply because, in case of the 1 -dimensional 
Hubbard model, analytical expressions for the holon (v/,) and spinon (v s ) velocities are given by 



M135H, 

v h = 2t sin (itn) ; v s = 



sin (27m) 
2%n 



(21) 



where, t and U are the nn hopping matrix element and the one-site Coulomb repulsion term, 
respectively, and n is the particle density (n < 1). Clearly, the velocity of the holon does not 
depend on U, while that of the spinon decreases, as established also from our tdDMRG studies. 
Furthermore, as U — > °°, spin velocity goes to zero. The holon moves by virtue of the hopping 

2 

matrix element while the effective spin-spin exchange, which is of the order of jj, propagates 
the spinon. And, in the thermodynamic limit, that is, U — > °° limit, only the holon propagates, 
the spinon doesn't "move" at all, v s being zero. As the magnitude of U increases, the velocity 
of the spinon decreases. 
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Time evolution of £t)> for Hubbard ctLainc Time evolution of *S CtJ* for Hubbard cbainc 



(max =200; u =2 .03: 6 - O.Or (doltsd curvs); 6 - 0.0 (solid curva) (max = 2.0O; u = Z.OJr 6- 0.0 r (dotlsd curve): 6- 0.0 (solid curva) 
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Figure 3: Time evolution of (Left) and (SfXO) (Right) for uniform (solid curve) and dimerized 

(dashed curve) Hubbard chains of length 10, 20, 30 and 40 sites, with = 2.0. Dimension of the DMEV 



basis used is 200. 
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'Dime evolution o£ "in^Ct)*- for Hubbard ctLainc 
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Time evolution of ^-S^ CtJ* for Hubbard cbaino 

Cmax = 2.0O; tt = i.oj r c — 0.0 7 (dotl sd curvs'l ; & — 0.0 (solid curva) 



0.5- 



j0.5- 



c 

V 



0.5- 



"V" 



4 0- site chain 



l i l I 



30-site chain 



20 -site chain 



i . i i ■ i . i ■ i 




i . i 



10- site chain 

i.ii,! 



o 

-0.2 
-0.4 


-0.2 
-0.4 



10 15 20 25 

time [fs) 



3 





-0.2 
-0.4 



-0.2 
-0.4 



40-sifce chain 



l.l.l 




2Q-3ite chain 



i ■ i i i 




10— aite chain 



■ ■■■■■■■■■ ■ 



A 



V 



10 15 20 25 

time (fs) 



30 



Figure 4: Time evolution of (Left) and (S^(t)) (Right) for uniform (solid curve) and dimerized 

(dashed curve) Hubbard chains of length 10, 20, 30 and 40 sites, with = 4.0. Dimension of the DMEV 



basis used is 200. 
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'Dime evolution of £fc]> for Hubbard clLainc 
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(dashed curve) Hubbard chains of length 10, 20, 30 and 40 sites, with = 6.0. Dimension of the DMEV 



basis used is 200. 
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Time evolution of <n L (t) > for Hubbard 
chains (max = 200; S = 0.0) 
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Figure 6: Comparison of the time evolutions of («l(0) (Left) and (S z L (t)) (Right) for uniform Hubbard 
chains of length 10, 20, 30 and 40 sites, with U, = 2.0 (solid curve), 4.0 (dashed curve) and 6.0 (dotted 
curve). Dimension of the DMEV basis used is 200. 
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Time evolution of <n L (t):> for Hubbard 
chains (max = 200 r 5 = 0.07) 
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Figure 7: Comparison of the time evolutions of (Left) and {S z L (t)) (Right) for dimerized Hub- 

bard chains (8 = 0.07) of length 10, 20, 30 and 40 sites, with ^ = 2.0 (solid curve), 4.0 (dashed curve) 
and 6.0 (dotted curve). Dimension of the DMEV basis used is 200. 



15 



5 Summary and Outlook 



To summarize, we find from our tdDMRG calculations that when a hole with a desired spin is 
injected at one end of the 7t-conjugated backbone of a polyene chain, it propagates from one end 
of the chain, to the other. The motion of this hole can be monitored by focusing our attention 
on the time evolution of the charge and spin densities at the two end-sites of the chain. In the 
absence of any external reservoirs (source-drain), the hole gets reflected back and forth across 
the length of the chain showing oscillatory motion. In the absence of electron-electron corre- 
lations, the charge and spin degrees of the hole do not get decoupled. The time taken by the 
hole to travel across the whole polyene backbone increases with approximately constant veloc- 
ity. We are currently extending these studies to PPP model and polymer topologies involving 
phenyl rings. It is seen that for dimerized chains, the velocity decreases further because of the 
fact that velocity of the hole is determined by the smaller of the two hopping matrix element, 
tii+l = to(l — (8)'), where 8 is the dimerization parameter and to is the mean hopping matrix 
element. For Hubbard chains, where spin-charge separation occurs, the hole "breaks-up" into 
two elementary excitations, one carrying only charge (holon), and the other, only spin (spinon), 
both of which moves with different velocities. It is found in accordance with the earlier litera- 
ture, the holon moves faster than the spinon, and with increasing U, although the velocity of the 
holon remains "almost" unaltered, that of the spinon significantly decreases. 
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